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Abstract. We consider a model of random permutations of the sites of the cubic 
lattice. Permutations are weighted so that sites are preferably sent onto neighbors. 
We present numerical evidence for the occurrence of a transition to a phase with 
infinite, macroscopic cycles. 
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1. Introduction 

Geometric representations of systems of statistical physics have a long history going back 
to the treatment of the Ising model by Peierls [§]. The Feynman-Kac formula provides 
such a representation for quantum models. It was originally introduced for the Bose 
gas jl], where the symmetric nature of the particles leads to the occurrence of random 
permutations. Siito showed in the ideal gas that Bose-Einstein condensation occurs if 
and only if infinite cycles are present |10[ ITT] . However, this relation does not seem to 
be always true in interacting systems [13]. Spin models also have representations where 
correlations are represented by loops, see e.g. [TJ [3] IT2] . 

The behavior of systems of interacting particles is notoriously difficult. On the other 
hand, the presence of infinite cycles should not depend closely on the microscopic details 
of the model. Results about simpler models of cycles could prove useful, especially for 
the understanding of the critical behavior. The purpose of the present article is to discuss 
such a model. 

An appropriate model of random permutations must involve the spatial nature of the 
original physical system. Namely, particles are spread over a certain domain, and distant 
particles are not directly correlated. The behavior of the model should also depend on the 
dimension of the space. Recall that the Bose-Einstein condensation takes place in three 
dimensions, but not in one or two. It seems therefore natural to consider permutations 
7L d — > 7L d instead of N — > N. For A C Z d finite, we will take the probability for a 
permutation tt : A — > A to be proportional to 

J~J e -a|x-7r(x)| 2 _ 

iceA 

The motivation for these Gaussian weights comes from the Wiener measure for bosonic 
trajectories in the Feynman-Kac representation of the Bose gas. The parameter a is 
proportional to the temperature of the system (not to the inverse temperature, as usual 
in statistical mechanics!). This model was progressively introduced in ® [5] . 

There should be no long cycles when a is large, i.e. when sites are heavily discouraged 
from jumping to a neighbor. Cycles should increase in size when a decreases. The main 
question is whether a transition occurs for some value a c > 0, below which a fraction of 
the sites find themselves in infinitely long cycles. 
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We present numerical evidence that points towards the occurrence of infinite cycles in 
three dimensions. A fraction of sites with positive density belong to infinite cycles. A 
natural question is how many large cycles are found in a typical permutation, on a cube 
of size LI One can argue that cycles are similar to a closed random walk, which has 
Hausdorff dimension two, so their number should grow like L 3 /L 2 = L. But one can also 
argue that infinite cycles represent correlations in the Bose condensate — two particles 
are correlated if they belong to the same cycle. In this case, there should be just one 
macroscopic cycle. 

The surprise is that neither conclusion above is correct: As it turns, infinite cycles are 
macroscopic, i.e. they involve a positive fraction of the sites, and their number fluctuates. 
This behavior was observed in the ideal Bose gas jl 1] . 

The model of random permutations is introduced in Section |2j where it is also shown 
that there are no infinite cycles at high temperature. A detailed probabilistic setting is 
described in Section |3j Our numerical results are presented in Section |4j They show a 
surprisingly close relationship between the random cycle model and the ideal Bose gas. 
The systems do not only share similar features qualitatively, but also quantitatively! 

2. The model 

As is usual in statistical mechanics, we define the model first in a bounded domain, and 
then consider an appropriate thermodynamic limit. Let A C Z d be a large but finite cubic 
box centered at the origin, and let Ba denote the set of permutations on A (i.e. bijections 
A — > A). The probability for tt £ Ba is given by 

^a(vt) — [J e-«l*-^)l 2 . (2.1) 
^ ' xeA 

Here, the parameter a represents the temperature, | • | denotes the Euclidean distance in 
and Z(A) is the partition function 

Za = Yl II e" a|a; - w(x)|2 . (2.2) 

7reB A x£A 

A cycle 7 of length k is a A;- tuple of distinct sites (x\, . . . , Xk). Given a permutation n, 
we let 70 denote the cycle that contains the origin, and £q = ^o( 7r ) its length. That is, 

70= [0M0)y(0),...y°-\0)). (2.3) 

Notice that 1 ^ £ ^ |A|. 

The interesting phenomenon is the possible occurrence of infinite cycles, that may oc- 
cur in the thermodynamic limit. This motivates us to introduce the function (p(a) that 
represents the probability that the origin belongs to an infinite cycle. First, we define 

7reB A ^o=fc 

and 

PU = k)= lim P A Uo = k). (2.5) 

Ayz d 

The existence of the thermodynamic limit is guaranteed by Cantor's diagonal argument: 
There exists a subsequence A m of increasing cubes such that PA m (^o = k) converges 
simultaneously for all k. In this paper, the notation A /* Z rf always refers to this specific 
subsequence. It is clear that ^2 k P\(£o = k) is equal to one for each finite A. But it 
may be strictly less than one in the limit A y Z d (th is is formalized by Fatou's lemma in 
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analysis). The probability (p(a) that the origin belongs to an infinite cycle is then defined 
by ' 

fc> i 

Let a c denote the critical temperature, 

q c = sup{a : ip(a) > 0}. (2-7) 

We expect that (p(a) is monotone decreasing in a (and strictly monotone decreasing for 
a < a c ), but we are unable to prove it. If one chooses a = in (2.1), one easily checks 
that P A (4 = k) = Aj for all k; then P(£ = k) = 0, and ^(0) = 1. 

It is easy to show the absence of infinite cycles when a is large enough. The following 
theorem implies that a c < oo for any dimension. 



Theorem 2.1. Assume that 



£ e-"W 2 < 1. 

x& d \{0} 



Then tp(a) = 0. 

Proof. Given a cycle 7 = (x\, . . . , Xk), let 



u>(i) = "[I e - a ^+ l ~ x ^ (2.8) 
j'=i 

(we set x/j+i = x\). From the definition (|2.1[), we have 



ft«. = *) = £ •l / IMI ^f"" P-9) 

7=(0,x 2 ,...,x fe ) 

The ratio of partition functions is less than one. Since 0^(7) ^ ITfci e~ a \ x j+ 1 ~ a: i\ 2 , and 
neglecting the restriction that cycles do not self- intersect, we get 

Pa(4 = *0 < ( E e- !'! 3 )*" 1 . (2.10) 
xez d \{o} 

Then P\{Iq = k) ^ c fe_1 with c < 1 for all domains A. From the definitions (2.4)-(2.6), 
we have 

A 

09(a) = lim lim 1 - V P A (£ = k) 

IM ^ 



lim lim V P A (£ = k). 



/ fc=A+l 



The sum is bounded by a convergent geometric series uniformly in A, see (2.10), which 
implies that (p(a) = 0. □ 

We can define a thermodynamic potential by 

/(a) = lim^logZA. (2.12) 



Here, we can take the limit along any sequence of boxes of increasing sizes, as can be estab- 
lished by a standard subadditive argument. The function /(a) is convex, since log Z\ 
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is given by the expectation of positive fluctuations. We conjecture that f(a) is analytic 
for all a, except for a = a c . 

3. The probability model 

This section is more technical and it can be skipped on first reading. But we hope 
to catch the interest of probabilists, which prompts us to turn the ideas above into a 
probability model. The situation is actually not trivial; our results require a condition, 
see (3.8) below, that we cannot prove in the case of Gaussian weights. The elements of 
the probability space are lattice permutations. The a- algebra contains the event that the 
origin belongs to an infinite cycle, and the probability measure is inspired by (2.1 ). 

3.1. The probability space. Let B be the set of all permutations ir on Z rf (i.e. bijections 
Z d -> Z d ). We redefine JB A so that it is a subset of B, by setting 

B A = {vr G B : vr(x) = x for all x <£ A}. (3.1) 
Let B xy be the set of permutations such that x is sent onto y, 

B xy = {ireM:Tr(x)=y}. (3.2) 
We let B denote the a-algebra generated by {B xy } xy£Z d. For any n and any x\, . . . ,x n , 
Vl, • ■ ■ j Vn S the probability of the set 



is defined by 



B = n B xm (3.3) 

x y ' vre-BnB A xeA 

Here, y) is a nonnegative, symmetric function on Z rf x Z d . An interesting example is 
i{x,y) = \x — y\ 2 . The normalization Z(A) is given by 

Z(A) = II e~ a ^ x ' n{x)) . (3.5) 
7reB A xeA 



The thermodynamic limit A / Z in (3.4) exists at least on a subsequence of increasing 
cubes, thanks again to Cantor's diagonal process. 

We now introduce the event that the origin belongs to an infinite cycle. First, consider 

Bi k) = U H B x . x ,, , (3.6) 

where x\ = x^+i = 0. The union is over distinct sites X2,...,Xf~ G Z rf (they are also 

(k) 

distinct from 0). The set B$ represents the event that the origin belongs to a cycle of 
length k. Next, let 



(DO) 



U B 

k > l u 



c 



(3.7) 



be the event for the origin to belong to a cycle of infinite length. It is clear that B^°^ 
belongs to the a- algebra B. 

The proof of the existence of a measure is not trivial, and we need the following as- 
sumption: For any x G Z d , we suppose that 
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This is equivalent to another condition that is more technical, but also more explicit: 

Pa(B w ) = 0. (3.9) 



lim lim 

\y-x\>R 

This condition means that sites do not jump straight to infinity in one step. Notice that 
it fails to be true when a = 0. It trivially holds if £(cc, y) = oo when \x — y\ is larger than 
some cutoff distance. It can also be established for £(x,y) = \x — y\ 2 when a is large, but 
we cannot prove it for arbitrary a > 0, although it is certainly true. 



Theorem 3.1. Assume that (3.8) holds true. Then there exists a unique probability 
measure p on (B,£>) that coincides with (3.4) for all sets of the form (3.3). In addition, 

we have p(Bq°°^) = ip(a), where (p(a) is defined as in Section 2, but with instead 
of\x-y\ 2 . 



This theorem is proved at the end of Section 3.2 The condition (3.8) is necessary. 
Otherwise, there exists x G Z d such that ^2 y p(B xy ) = 1 — e with e > 0. Consider then 

Then n n B n = 0, but 



the decreasing sequence of sets (B n ) with B n = 

p(B n ) = 1 - Yl 



^y,\y—x\>n^xy 



p{B xy ) ^ e, 
y,\y-n\>n 

so p(B n ) does not go to and p cannot be a probability measure. 



(3.10) 



Notice that Theorem 2.1 in Section [2] extends straightforwardly to general £(x,y). It 

,(o 

o 



follows that P( J B l (oo) 



when a satisfies Ylx^=o e a ^°' x ' < 1. 

3.2. Construction of the measure. Let us introduce the algebra B' generated by the 
sets B xy . It is not hard to verify that the set function defined in (3.4) extends uniquely 
to a finitely additive measure on B' . We need to prove that it is a- additive within the 
algebra. 



Lemma 3.2. Suppose that the property (|3.8|) is satisfied. Let B\ D By, D . . . be any 

Pill, 



Then lim r 



0. 



sequence of sets of B' that decreases to 0, i.e. n n ^i B n 

Proof. We show the counterpositive, namely that if (B n ) is a decreasing sequence in B' 
such that p(B n ) > e > for any n, then r\ n B n is not empty. 

Given n G B and A C Z rf , we denote A(tt] A) the set of permutations whose restriction 
to A coincides with it. Precisely, 

A(ir; A) = {vr' G B : vr'(x) = vr(x) for all x G A}. (3.11) 

Since B n G B' , for each n there exists a finite set A n such that B n is given by 

B n 



U A{ir;A % 



(3.12) 



The set B n is uncountable, but there are only countably many distinct sets A(ir; A n ). Thus 
Eq. (3.12) is really a countable union of disjoint sets. We can suppose that (A n ) is an 
increasing sequence of cubes centered at the origin. Let k n be some integer, and C n C B n 
be the set of permutations where each site of A n is sent at distance less than k n : 

C n = {tt £ B n : |vr(x) — x\ ^ k n for all x G A n }. (3.13) 

We have 



B n \C n c U 



U B 



xy 



(3.14) 



\y— x\>k„ 
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It follows that 



p(B n \ C n ) = lim p A (B n \ C n ) 



< y lim y Pk{B xy ). 



x6A r , 



(3.15) 



From the assumption (3.9) we can choose k n large enoug h so that p(B n \ C n ) < |e. Then 
(C n ) is a decreasing sequence such that p(C n ) ^ \e for all n. 



In a fashion similar to (3.12), we can decompose C n into the disjoint union 



C n = U A(TT n<i ;A n ) 
i=l 



(3.16) 



for some suitable permutations 7r nj j. The number r n is now finite. We have A n C A n +j. 
Two sets A(ir,A n ) and j4(vt', A n +i) are either disjoint, or A(tt, A n ) D A(7r', A^+i) if the 
restrictions of 7r and ir' on A n coincide. We can define a tree with vertices (re, i), n = 
1,2, .. . and 1 ^ i ^ r n , and with an edge between (re, i) and (re + l,j) whenever 



A(lT nt i,A n ) D A(7T„ + 1J, A n+ i) 



(3.17) 



We also connect the vertices 1 ^ i ^ ri, to the root 0. There are r n vertices 

at distance n from the root, so all incidence numbers are finite. For each vertex (n, i), 
let q(n,i) denote the length of the longest path descending from (n, i). For each n it is 
infinite for at least one i. In addition, each (re, i) with q(n, i) = oo is connected to at least 
one vertex (re + l,j) with q(n + 1, j) = oo. These properties hold true because incidence 
numbers are finite, i.e. because we defined the tree starting with the sequence (C n ) instead 
of (B n ). We can select an infinite path (n,j n ) with q(n,j n ) = oo for all n. Since any x 
belongs to some A n , we can define a map 7r : T, d — > Z rf by setting 



(3.18) 



Observing that ir is a permutation, and that 7r G C n for all re, we see that D n C n is not 
empty. □ 



Proof of Theorem 3.1 It follows from Lemma 3.2 that p is a cr-additive premeasure on 
B'. It has a unique extension to a measure on ,6 by the Car atheo dory- Frechet extension 
theorem. 

In order to show that p(B^°^) is equal to the probability of infinite cycles as defined 



(k) 

in (2.6), let us observe that B\ fl B { 



if k + k', and B^° ] = [U fc S^ ; ] c . One can 

Q fc ^). The result then follows from the 

□ 



check that P(£q = k) defined in (2.5) is equal to p(B t 
cr-additivity of p. 



4. Numerical results 

4.1. Description of the method. We have performed intensive Monte Carlo simulations 
of the random cycle model with Gaussian weights for the jumps. The dynamics is pure 
Metropolis; a change in the permutation configuration is accepted or rejected according 
to the change in "energy". 

A "step" of our code consists in sweeping the sites in lexicographic order. For each 
site x, we randomly choose a site y in a window centered at x with size depending on the 
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temperature a. Given x, y, we consider changing the permutation tt into tt', where ir 1 is 
defined as follows (see Fig. [T]) : 



vr'(x) = y, 



(4.1) 



and tt'(z) = ir(z) for z ^ x, tt (y). 



X 



7T[X , 



y 



Figure 1 . Illustration for the update scheme of the algorithm. Bold lines 
represent the old permutation tt, dashed lines represent the new permuta- 
tion 7r'. 

The energy difference between new and old permutations is AH = H\{tt') — H\(tt), 
where H\{tt) = X^xeA \ x ~ ^C 3 -)! 2 denotes the "Hamiltonian" of the model. We have 

I i2i_i i2i i2i_i i2 

AH = \x — y\ + | vr (y) — vr(x)| — \x — 7r(x)\ — | vr (y) — y\ . (4.2) 
Then according to the Metropolis prescription, the change is accepted with probability 

min[l, e~ aAH }. 

The initial configuration is usually the identity permutation ir(x) = x, but we have also 
considered initial random configurations chosen over all permutations A — > A with uniform 
distribution. It turns out that the system thermalizes extremely well, irrespective of the 
initial condition. Measurements are taken after suitable thermalization. All numerical 
computations were performed on a personal computer and they can be reproduced rather 
easily. 

Most of our measurements are for the random variable pj-, that represents the fraction 
of sites that belong to cycles of length less than or equal to k. Precisely, pk is defined by 

#{x G A : 4(vr) ^ k} 



|A| 



(4.3) 



where £ x is the length of the cycle that contains x. We have ^ Pki^) ^ 1 and 1 ^ k ^ |A|. 
It is related to </?(«) if we assume that the system is essentially translation invariant. We 
expect that 

¥>(<*) « 1 - (p N ), (4.4) 

with N such that 1< N < |A|. 
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4.2. Fixed temperature, different dimensions. We have first fixed a = 0.2 and con- 
sidered cubic boxes in dimensions d = 1, d = 2, and d = 3. Fig. [2] depicts the graphs of 
the expectation {p^} of the fraction of sites in cycles of length less than k. The horizontal 
axis in Fig. [2] are fe/|A|; they take values between and 1. 





0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



(a) 



(b) 



(c) 



Figure 2. Expectation value for p^ as function of k/\A\ for a = 0.2. (a) 
d = 1 and |A| = lO'OOO. (b) d = 2 and |A| = 100 2 . (c) d = 3 and |A| = 50 3 . 

Fig. [2] (a) shows that almost all sites belong to cycles with very small length compared 
to the volume of the system. The situation is different in d = 2 and d = 3. In d = 2, 
around 25% of the sites belong to small cycles, and 75% belong to macroscopic cycles, i.e. 
to cycles whose length is a fraction of the volume. We expect that this density is equal to 

3, with respective densities 3% and 97%. 

2, where there are clearly no 



(f(a) defined in (2.6). The same holds for d 
Fig. [3] shows the situation at the higher temperature a 
macroscopic cycles. 



0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



(a) 



(b) 



(c) 



Figure 3. Expectation value for p^ as function of k/\A\ for a = 2. (a) 
d = 1 and |A| = lO'OOO. (b) d = 2 and |A| = 100 2 . (c) d = 3 and |A| = 50 3 . 



From these preliminary exploration, it seems that macroscopic cycles are present in 
dimensions greater or equal to 2, instead of 3! 

4.3. Two dimensions. The main numerical result for d = 2 is shown in Fig.|4j where the 
expectation (p^.) is plotted for a = 0.1 and different lattice sizes. It is manifest that finite 
size effects are very important. The density of sites in small cycles grows from 10% for 
L = 50 to 30% for L = 2000. We expect the curves to continue their progression upwards 
as the size increases, until no macroscopic cycles is left in the limit L — > 00. 

A possible explanation involves random walks. The cycle containing the origin is a 
self-avoiding closed random walk — but its probability differs from that of random walks 
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fc/|A| 



Figure 4. Expectation (pk) for d = 2, a = 0.1, and with various domain sizes. 



because of the presence of all other cycles. Nonetheless, the analogy with random walks is 
worth pursuing. Random walks are recurrent in d = 2. The probability f n for the simple 
random walk to return to the origin for the first time after n steps satisfies j2] 



n > 1 



fn 



1. 



f n ~ (nlog 2 n) J 



(4.5) 



The "macroscopic cycles" in d = 2 are those that are big with respect to the size of the 
domain, i.e. whose mean square distance y/n is larger than the size L of the domain. Let 
p(L) denote the density of sites in long cycles. The condition y/n ~ L implies that the 
probability for the origin to belong to a long cycle is roughly 

-i 



P(L) 



i > L' 2 



(logL)- 



It follows that 



p{L)-p{2 k L) 



(4.6) 



(4.7) 



In first approximation this is proportional to k, as observed in Fig. |4j 

To summarize this section about two dimensions, we have understood that there are no 
really macroscopic cycles, and that the evidence suggested in Fig.[2](b) comes from strong 
finite-size effects. 

4.4. Three dimensions. Fig. [5] shows (pk) for a = 0.8 and different sizes of the do- 
main. In contrast to Fig. [4] there are no noticeable finite-size effects. It is thus clear that 
macroscopic cycles are present in three dimensions, and that numerical simulations work 
remarkably well. 

We have computed the density of sites in macroscopic cycles as a function of the tem- 
perature, see Fig. [6] Recall that it should be equal to (f(a). We find that ip(a) seems to 
be continuous, and that the critical temperature is a ~ 1.7. It is instructive to compare 
it with the critical temperature for the Bose-Einstein condensation of the ideal gas, which 
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is known exactly. The difference between our model and the Feynman-Kac representation 
of the ideal gas is that our particles are frozen on the sites of the cubic lattice. But the 
comparison is otherwise possible. The particle density in our system is equal to one, and 




Figure 6. Density of sites in macroscopic cycles versus temperature, in 
three dimensions. The points represent numerical data for L = 50 (in 
black) and L = 100 (in white). We observe that the critical temperature of 
our model is a c ~ 1.7. The curve gives the density of the Bose condensate 
p$ EC (a), Eq. p|. 



the physical constants are ^ = k-g = 1. Then our a is proportional to the temperature, 
with a = 1/4/3. With £ denoting the Riemann zeta function, the critical temperature of 
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the corresponding ideal Bose gas is 

c^ EC = C(3/2r 2/3 7r = 1.656... (4.8) 

It is tantalizingly close to the critical temperature in the random cycle model! In addition, 
the density of the condensate in the Bose gas is given by 

^ EC («) = 1-(^ec) 3/2 - M 

It is plotted in Fig. [6] along ip(a). (f(a) and p® EC (a) appear to be close, but it seems that 
their differences cannot be accounted for by numerical errors of by finite size effects. 

All numerical results so far were for the average (pk) over many permutations. We now 
discuss the distribution of macroscopic cycles in a typical permutation. Fig. [7] shows pk 




27000 



1 10000 27000 1 10000 27000 

Figure 7. Thin lines show pk{^) in the permutation after 1, 100, 1000, 
and 20'000 steps. The system has been thermalized initially. Bold lines 
show the average {pk) in all permutations until 1, 100, 1000, and 20'000 
steps. Here, d = 3, a = 0.5, and L = 30. 

for a = 0.5 and L = 30. We plot the value of pk{^) for the permutation tt after r steps, 
and its average {pk) over all permutations before r steps. The latter converges as r — > oo, 
as expected. The jumps in the graph of pk(n) correspond to macroscopic cycles. A jump 
at tttt means that a macroscopic cycle of density m is present, and accordingly the height 
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Q 


/ ^max \ 

\ |A| / 


ip{a) 


<lAf>M«) 


0.001 


0.6242 


0.9998 


0.6243 


0.25 


0.5804 


0.9295 


0.6244 


0.50 


0.5115 


0.8187 


0.6247 


0.75 


0.4338 


0.6944 


0.6248 


1.50 


0.1414 


0.2264 


0.6244 


1.68 


0.0251 


0.0402 


0.6258 



Table 1. Average density of the longest cycle as a function of the tem- 
perature. Here, d = 3 and L = 50. 



of the jump is given by 

p k (ir) - Pk-^n) = fa. (4.10) 

The number of macroscopic cycles (of density larger than e > 0) seems to fluctuate. As 
A y Z 3 , the distribution of this number converges to some nontrivial distribution. This 
behavior was observed by Siito in the ideal gas [11] . Macroscopic cycles are due to particles 
in the condensate, for which the underlying distribution of permutations is uniform. Let 
us tentatively extrapolate this observation to our model. It suggests that the domain 
A splits into two sets Af and A^, that corresponds to sites in finite and macroscopic 
cycles, respectively. While these sets are random, their respective volumes are close to 
(1 — (^(a))|A| and ^(ck) |A|. Suppose that the distribution of macroscopic cycles in Aqo is 
the same as if the permutation on A^ was chosen with uniform probability. The average 
density of the longest cycle in uniformly distributed random permutations is known, see e.g. 
j9], and is equal to 0.6243. . . . We report in Table [T] the results of numerical measurements 
in our model. We find values that seem to be in total agreement. This is extremely 
surprising, since the jumps in our permutations satisfy certain spatial restrictions; with 
uniform permutations, sites can go to infinity in but one step. 



5. Conclusion 

We have considered a model of random permutations on the cubic lattice. It is a 
priori a crude approximation for the ideal Bose gas in the Feynman-Kac representation. 
Surprisingly, its behavior is close to the Bose gas both qualitatively and quantitatively. The 
critical temperatures are very close, and the average densities of the longest cycle coincide. 
But there seems to be a small difference between the density of sites in macroscopic cycles 
and the density of the Bose-Einstein condensate. The study of larger systems on powerful 
computers may shed more light on this issue. 

The fact that this simple model is so close to the ideal Bose gas suggests to use it in 
order to simulate interacting systems. A natural generalization is to introduce interactions 
between permutation jumps. Can we do it so that the critical temperature of the random 
cycle model is in quantitative agreement with the critical temperature of the interacting 
Bose gas? 

Many mathematical aspects need clarifying as well. One would like to know about 



analytic properties of the thermodynamic potential (2.12). The existence of a probability 
measure on the space of permutations on Z d should be completed. The absence of cycles 
in one and two dimensions should be established, and also certain properties such as the 
monotone decreasing behavior of ^p(a) with respect to a. 
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